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The alignment of molecular levels with the Fermi energy in single molecule junctions is a crucial 
factor in determining their conductance or the observability of quantum interference effects. In the 
present study which is based on density functional theory calculations, we explore the zero-bias 
charge transfer and level alignment for nitro-bipyridyl-phenyl adsorbed between two gold surfaces 
which we find to vary significantly with the molecular conformation. The net charge transfer is the 
result of two opposing effects, namely Pauli repulsion at the interface between the molecule and the 
leads, and the electron accepting nature of the NO 2 group, where only the latter which we analyze in 
terms of the electronegativity of the isolated molecules depends on the two intra- molecular torsion 
angles. We provide evidence that the conformation dependence of the alignment of molecular levels 
and peaks in the transmission function can indeed be understood in terms of charge transfer for this 
system, and that other properties such as molecular dipoles do not play a significant role. Our study 
is relevant for device design in molecular electronics where nitrobenzene appears as a component in 
proposals for rectification, quantum interference or chemical gating. 

PACS numbers: 73.63.Rt, 73.20.Hb, 73.40.Gk 



I. INTRODUCTION 

In the emerging field of single molecule electronics^ 
considerable progress has been made in recent years 
in achieving comparable results for the electron trans- 
port through nano junctions from theory and experi- 
ment for a small range of well-defined test systems such 
as H2 molecules between Ft leads^-^ or bipyridin o ^ — P 
and alkane mono- or di-thiolates^^-^^ between Au elec- 
trodes. None of these systems, however, show interesting 
device properties and although a variety of theoretically 
designed single molecule devices^^-^^ exist in the litera- 
ture, these proposals usually assume an idealized atomic 
configuration for the junction setup, where most degrees 
of freedom for the position of the nuclei are frozen and 
which can rarely be matched by experiments. A molec- 
ular component, which is particularly popular for molec- 
ular rectifiers^l^ but also more recently has been sug- 
gested in the context of quantum interference effect o^^ 
is a conjugated or cross-conjugated tt system substitued 
with one or more NO2 groups. 

It is by now well established that the alignment of the 
eigenenergies of molecular orbitals (MOs) with the Fermi 
energy of the metal electrodes^^-^^, which is a key fac- 
tor in determining the conductance of a molecular nano 
junction, is strongly related to the zero-bias charge trans- 
fer between the molecule and the leads^^-^^ where also 
the effect of Fauli repulsion at the interface (often ad- 
dressed as interface dipole or pillow effect) can play a 
^Qlg33,34 Hence, the electron accepting nature of nitro- 
substituents is routinely used to induce an upwards shift 
of molecular levels^^ which is also supported by theoret- 
ical trend studies for ideal configurational setups of the 
junction geometry^^^^^. Such shifts have even been used 
to propose the concept of chemical gating^^, where data 
has been presented for a variety of electron accepting and 



donating groups but NO2 has not been included in the 
study. In a recent work based on photoemission spec- 
troscopy^, however, empirical evidence was found, that 
the assumed strong correlation between the donor or ac- 
ceptor strength of chemical substituents on the one side, 
and level alignment or work function reduction (in the 
case of adsorbed monolayers) on the other side does not 
always hold, and hints at a dependence of this correlation 
on molecular orientation were given^. While the depen- 
dence of the electronic coupling on molecular conforma- 
tion in biphenyl with rather strong coupling to gold leads 
has been investigated recently both by experiment^^^^^' 
and theory^^^^^'^^, no attention has been paid so far to 
the effect of finite torsion on level alignment, which is 
more decisive for the conductance in systems with weaker 
bonds to the electrodes as it was found in a comparison 
between bipyridine and biphenyl-dithiolat o^^i^^ . 

In the present work we show in calculations based 
on density functional theory (DFT) and employing the 
Siesta code^^ that for nitro benzene anchored to Au elec- 
trodes by pyridyl linkers (Fig. [TJi) (a system we recently 
investigated in the context of quantum interference ef- 
fects^^) the molecule loses fractions of electrons to the 
leads for a wide range of geometries in spite of the elec- 
tron accepting nature of the NO2 group. This finding is 
explained by i) relatively low values for the electron affin- 
ity (EA) of nitro-bpph and ii) the influence of Fauli re- 
pulsion effects at its interface to the leads as a mechanism 
for electron removal from the molecule which we explored 
in great detail in earlier work for the Au/bipyridine/Au 
system^^^^^. Although the values for EA are rather low, 
their conformation dependence still determines the de- 
pendence of charge transfer on the two torsion angles 
and (j)2 in Fig. [1^, where we have studied 100 differ- 
ent atomic configurations for the molecule in the junc- 
tion, while Fauli repulsion at the interface is found to 
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be conformation independent. We also present evidence 
suggesting that in this junction setup the dependence of 
the molecular level alignment on and 02 is dominated 
by the variation of charge transfer (which includes in- 
terfacial Pauli repulsion in our definition) and that for 
this dependence the molecular dipole and changes in the 
composition of MOs should play a rather insignificant 
role. 

The paper is organized as follows: In the next section 
we describe our computational methods and introduce 
the scope of molecular conformations studied, where we 
also address the question of their respective stability. The 
following section contains a discussion of our results in 
five steps. First we establish the conformation depen- 
dence of transmission functions and molecular eigenen- 
ergies in Au/nitro-bpph/Au junctions which is the main 
motivation of our study and also its link to our previous 
work in Ref.^^. Secondly, we discuss equilibrium charge 
transfer and its driving forces for the single molecule 
junctions in our article, where the effects of both interfa- 
cial Pauli repulsion and molecular electronegativity are 
discussed. In a third step we argue that charge transfer 
explains the conformation dependence of the alignment 
of molecular levels rather than molecular dipoles or a 
change of composition in MOs. We then check explic- 
itly the impact of shortcomings of our DFT approach on 
the conclusions in this article by applying ad hoc cor- 
rections for the underestimation of energy gaps and the 
methods failure to describe dynamical screening effects. 
Finally, we describe the changes in molecular dipole mo- 
ments with the two torsion angles 0i and 02- The last 
section contains a summary of our main conclusions. 



II. METHODS AND STRUCTURES 
A. Computational details 

Throughout this article we have generated our data 
by using the density functional theory (DFT) based 
Siesta code^^ with a double-zeta polarized (DZP) ba- 
sis set, Troullier-Martins type norm-conserving pseu- 
dopotential^ and a Perdew-Burke-Ernzerhof (PBE) 
parametrization^ for the exchange correlation (XC) 
functional. For the electron transport problem a non- 
equilibrium Green's function (NEGF) approach^^ has 
been employed which in the details of its combination 
with Siesta has been described in Ref."^^. In order to ap- 
ply this method it is necessary to divide single molecule 
junctions into three regions, namely a left lead, a right 
lead and a scattering region. For all three regions in- 
dependent supercell calculations with periodic boundary 
conditions are performed, where in the plane perpendicu- 
lar to the transport direction we chose 3x3 sections of Au 
fee (111) to form the unit cell for all geometries in this ar- 
ticle. In a second step the leads are added to the Green's 
function of the scattering region as self energies allowing 
for the calculation of transmission functions within the 



framework of the Keldysh formalism^^. For the k-point 
sampling in the transverse plane a 4x4 grid has been used 
for all systems but it has been checked for one structure 
that 6x6 and 8x8 grids do not change our results in any 
way. 

A possible source of concern might be seen in our 
LCAO basis set which was shown recently to give inac- 
curate results for various properties of noble metal sur- 
faces such as surface energies and work functions^^. The 
reason for this deficiency was identified in the too rapid 
decay of the surface wavefunctions into the vacuum. This 
problem results only in slight errors of ^ 0.1 eV for sur- 
face energies but nevertheless can be dramatic in terms 
of percentages because the latter quantity is generally 
small for noble metals. The work function is found to be 
underestimated by 1 eV by a DZP basis set for gold 
in Refj^, since it is calculated directly from the electro- 
static potential in the vacuum which naturally depends 
very strongly on the quality of the description of the long- 
range decay of the wavefunction. Most importantly for 
our study, the DZP-LCAO wavefunctions were only seen 
to deviate substantially from accurate plane wave cal- 
culations beyond a distance of 2-3 A . In our work half 
the bonding distance between the molecule and the leads 
must be seen as the relevant length scale for wavefunction 
decay regarding e.g. the accuracy of interface dipoles. 
This length is 1.06 A and we therefore expect our main 
results to be unimpaired by the named difficulties. We 
concede, however, that these technical problems affect 
our description of the free gold surface and the decay of 
its electron density into the vacuum, which can have an 
influence on the absolute values for the transfer of frac- 
tional charges which we discuss in Sec. IIIIBI But such 
an error is clearly independent of molecular conforma- 
tion since the free Au surface is always the same for all 
systems in this article and therefore has no relevance for 
our discussion. We also stress that we use trends in frac- 
tional charges only as a means of elucidating the physical 
driving forces for level alignment and that such fractions 
of charges are not themselves observable quantities which 
could be measured in experiments directly. 



B. Molecular junctions and their stability 

Fig. [TJi shows the geometry of 2,5-bis(4- 
pyridyl) nitrobenzene coupled to Au (111) surfaces 
with an ad-atom where the N-Au distance on both 
sides is 2.12 A which has been identifled as equilibrium 
distance for this system in earlier work^^^^. In Refs.^^i^S 
the choice of this particular surface structure with 
an ad-atom has been motivated as due to its effect 
of bringing the lowest unoccupied MO (LUMO) close 
to the Fermi energy in the junction where a detailed 
comparison with flat surfaces has also been undertaken. 
Throughout this article the abbreviation nitro-bipyridyl- 
phenyl (nitro-bpph) has been used for the molecule for 
the sake of simplicity. The two torsion angles and 
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torsion angle (|)i(°) 

FIG. 1: (Color online) a) Geometry of the Au/nitro-bpph/Au 
junctions in this article, b) Total energies for all angles 01 and 
02 with respect to the most stable conformation (01=40°, 
02=50°). The torsion dependence of the distance do-H be- 
tween the H and O atoms causing the largest contribution to 
steric repulsion (which are marked in a) ) is also displayed as 
inset. The green lines give an indication which structures are 
stable or could be stabilized by chemical means, namely the 
ones with E*'^* less than 1.5 eV or do-H above 1.2 A . 



02 have been varied in our calculations resulting in 
a total of 100 different molecular confomations. The 
atomic configuration of the molecule has been optimized 
once by DFT total energy minimization and then the 
angles have been varied rigidly. For the Au surfaces 
the atomic positions have been based on truncated bulk 
values. No relaxations of nuclei positions resulting from 
the interaction between molecule and leads have been 
considered in this article, because as can be argued from 
the rather small bonding energies for pyridil anchors 
on gold surfaces^, such relaxations could only result in 
minor corrections but this would come at the expense of 
increasing the computational effort substantially. 

A main purpose of our article is to investigate the in- 
fluence of the torsion angles 0i and 02 in nitro-bpph (see 
Fig. [TJi) on charge transfer and level alignment. In this 
section we ask the question how this two degrees of free- 
dom affect the stability of the molecule in the junction. 
For unsubstituted molecules consisting of directly con- 
nected aromatic rings such as biphenyl^^ and bipyridine^ 
it is well known that their tilt angles, ~ 45° and ~ 25°, 
respectively, are a result of the interplay between the tt 
conjugation, which tries to flatten the structure, and the 



steric repulsion of the hydrogen atoms in the ortho posi- 
tion with respect to the ring-connecting carbons, which 
favors a perpendicular arrangement of the planes of two 
adjacent rings. By chemical means this balance has been 
shown to be adjustable to lower angles by introducing 
a -CH2- bridge and to higher angles with bulky side- 
groups ^. 

For nitro-bpph the situation is more complicated, 
partly because there are now two conformational degrees 
of freedom (0i and 02) and also because the close prox- 
imity of one of the oxygen atoms of NO2 to the near- 
est hydrogen atom on the neighbouring pyridil results in 
an enormous energy barrier for planar or nearly planar 
conformations. The effect of this large energy barrier is 
demonstrated in Fig. [T]3 where we show that the total en- 
ergies of the eight conformations closest to planarity are 
1.5-15.0 eV higher than the one of the most stable confor- 
mation (01=40°, 02=50°). This result is not surprising 
when the values of do-H in the inset of Fig. [TJd are in- 
spected, where for these geometries the O and H atoms 
in question are found to be less than 1.2 A apart in dis- 
tance. Out of the remaining 92 conformations we studied 
in this article, however, an energy barrier of less than 0.5 
eV was calculated for 85 geometries, which puts them 
well within the range where chemical flxation techniques 
as demonstrated in Refs.^^^^^ can be employed for 0i or 
where the conformation might be influenced by an exter- 
nal electric fleld^^, a method that would be applicable in 
principle for both, 0i and 02- Without any form of flx- 
ation or external influence, both angles would vary over 
a wide range even for rather low temperatures in actual 
experiments due to the rather shallow energetic minima 
found in Fig. [TJ). We stress that the trends we derive 
and discuss for the conformation dependence of charge 
transfer and level alignment in our article are focused 
on the 92 structures that are experimentally achievable 
and results for the geometries close to planarity are only 
provided for the sake of completeness. 



III. RESULTS AND DISCUSSION 

A. Conformation dependence of the transmission 
function 

In Fig. [2^ we show the transmission function for two 
selected molecular conformations in order to illustrate 
that a variation of 0i can shift its peak structure close to 
the Fermi energy Ef signiflcantly thereby having a big 
impact on the conductance of the junction. We already 
discussed the consequences of the variation of 0i for the 
observability of quantum interference effects for this sys- 
tem in an earlier work^^, while the focus in the current 
article is on exploring the reasons behind the peak shifts 
which are a reflection of the torsion induced shifts of the 
MO onsite energies for the LUMO and LUMO+1 as dis- 
played in Fig. [JJd. For nitro-bpph there are in principle 
several explanations thinkable for the torsion angle de- 
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FIG. 2: (Color online) a) Transmission functions at energies close to Ef for two geometries with (01=0° 702=0°, solid line) 
and (0i=4O°/(/)2=O°, dashed line), b) Evolution of MO eigenenergies with 01 for the LUMO and LUMO+1 for 02=0°. The 
spatial distributions for these two orbitals are shown as insets for (0i=O°/02=O°). 



pendence of MO onsite energies: i) A change in confor- 
mation is a substantial variation of the molecular struc- 
ture, and therefore has an effect on the way MOs are 
formed as linear combination of atomic orbitals (AOs), 
which in turn might alter both the energetic distances be- 
tween different MOs and their spatial distribution, where 
the latter can have an impact on the couphng to the sur- 
face^ i^^i^^ . ii) Since the molecule is asymmetric, it ex- 
hibits a dipole moment which again might change with 0i 
or 02 and can play a role in the vacuum level alignment 
between molecule and leads, iii) As found for bipyridine 
between Au surfaces^^^^^, zero-bias charge transfer is ex- 
pected to be a significant factor for level alignment also 
for nitro-bpph and studying its variation with the molec- 
ular conformation is therefore an important issue. It is 
not straightforward to separate the above named effects 
in a rigorous way and in the following sections we will at- 
tempt such a separation and argue that charge transfer is 
the dominant factor for the variation of the MO energies 
with 01 and 02 in the Au/nitro-bpph/Au system where 
we base this conclusion on trends in quantities we can 
derive directly from our DFT calculations. 



B. Zero-bias charge transfer and electronegativity 

The starting point for our analysis has to be the ques- 
tion how many fractions of electrons are exchanged be- 
tween the molecule and the gold surfaces at zero bias 



for establishing an equlibrium in terms of chemical po- 
tentials. In Ref.^^ we addressed general aspects of this 
question. Fractional charges are only physically mean- 
ingful in the sense of describing the charge accumulating 
in the immediate vicinity of a molecule adsorbed between 
two metal surfaces, and do not necessarily describe the 
filling or emptying of MOs, because these MOs are hy- 
bridizing with the metal surface states and the system 
cannot be partitioned into its components in an unam- 
biguous way. The projection scheme we use in this article 
to derive energies for the LUMO and LUMO+1 is a way 
of performing such a partitioning but aimed at a special 
purpose. It is based on decoupling AOs on the molecule 
from AOs on the metal leads in the NEGF-DFT formal- 
ism, which is useful for identifing the molecular levels and 
their energies corresponding to peaks in the transmission 
function but a priori provides no information about the 
local charging of the molecule. 

Such information (as we established in Ref.^^) can only 
be reliably retrieved in two ways: i) The lowest molecular 
orbital (we call it MOl in the following) lies more than 
10 eV below all gold states and therefore is a pure molec- 
ular state, meaning that there is no hybridisation with 
the leads. The shift in its energy due to the interaction 
with the surface can therefore be only due to charging 
of the molecule, where the level goes up in energy when 
fractions of electrons are added and down if they are 
subtracted. A comparison of the MOl positions in the 
junction and for the isolated exphcitly charged molecule 
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then allows to calculate the local charges in an indirect 
fashion, ii) For rather weak covalent bonds (such as the 
ones formed with gold by pyridil anchors as was explic- 
itly demonstrated for Au/bipyridine/Au in Ref.^^), the 
integral of the charge density differences over the spatial 
region of the molecule AN^^f , where its density and the 
one of the Au slab are subtracted from the density of the 
combined system, also describes the charging very accu- 
rately because the border between the two systems can 
be drawn rather easily where the electron density has a 
minimum. 

We stress that charging as we define it here includes 
the effects of Pauli repulsion at the interface which are 
sometimes described as interface dipoles. Since the lat- 
ter mechanism pushes fractions of electrons away from 
the molecule, it results in effective electron removal and 
pushes down the energy of all MOs in the same way 
as emptying the highest occupied MO (HOMO) would. 
This was shown in detail for Au/bipyridine/Au in Ref.^^ 
where Pauli repulsion was identified as dominant mecha- 
nism for the level alignment. It was also shown explicitly 
that even the shape of the charge density differences is 
the same in emptying the HOMO for the isolated bipyri- 
dine and Pauli repulsion for the Au/bipyridine/Au inter- 
face. For Au/nitro-bpph/Au the same correspondence 
can be found, because also for this molecule the HOMO 
is localized predominantly on the pyridil anchors. 

In the top part of Fig. [3] we plot AN^^f in its full de- 
pendence on and 02 for the optimal bonding distance 
do=2.12 A . In spite of containing NO2 as a strong elec- 
tron acceptor, the molecule loses fractions of electrons to 
the surface for a wide range of angles. The latter net 
result is a consequence of the interfacial Pauli repulsion 
discussed above. This can be illustrated by comparing 
the distance dependence of AN^^f for (^i =407^2=0°) 
(dashed fine in the inset of the top part of Fig. [3]) with 
the study on Au/bipyridine/Au in Ref.^^. In both cases 
the fractional number of electrons on the molecule goes 
down when its distance to the electrodes is reduced from 
large values to do, because the repulsion effect decays 
continuously with an increase of the Au-N distance. A 
comparison of ANl^f for do (as shown in Fig. ^ and 
do +3 A (not shown here) reveals that the latter effect 
is rather conformation independent at least for values 
of and 02 higher than 30°. At do +3 A the lowest 
value for AN^^f (for both angles higher than 10°) is 
at (01 =90° /02 =90°), where it becomes zero. We also 
want to point out that the net charges on the nitro-bpph 
molecules are by a factor of ~ 4-5 smaller than those for 
the unsubstituted bpph (not shown here) or bipyridine 
as analyzed in Ref^^ which shows that there is a bal- 
ance between the electron accepting action of NO2 and 
interfacial Pauli repulsion explaining the relatively small 
amount of charging in the Au/nitro-bpph/Au system. 

For a further interpretation of the trends in AN^^f ^ 
it is useful to look at the torsion dependence of a re- 
lated property for nitro-bpph, namely its electronega- 
tivity EN. According to Mulliken^ it can be defined as 
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FIG. 3: (Color online) (Top) Variation of the transfer of elec- 
trons from the Au surfaces to the molecule ANl!^f derived 
from integration over charge density distribution differences 
between the whole junction and the molecule and surface in 
isolation^^ at bonding distance do (2.12 A ) in dependence on 
the angles 01 and 02, where the x-axis denotes 01 and differ- 
ent values for 02 are distinguished by varying line types. The 
inset shows the distance dependence of AN^^^ for the two 
angle combinations also highlighted in the transmission func- 
tions of Fig. [2^. (Bottom) Variation of the electronegativity 
EN of the molecule in dependence on 0i and 02 with its two 
components, EA and IP, displayed in the insets. 



EN=(IP+EA)/2, where the ionisation potential (IP) and 
EA can be calculated from DFT total energies of the 
charged and neutral molecule as IP=E(N-1)-E(N) and 
EA=E(N)-E(NH-1), respectively^. The bottom part of 
Fig. [3] reveals that at least qualitatively the conformation 
dependent trends in EN match those in AN^^f . From 
the insets it can be seen that the abrupt downward trends 
for low angles are mostly due to IP, while EA exhibits a 
continuous increase when going from high to low angles 
over the whole conformational range. 

The latter finding would also be expected from chemi- 
cal intuition, since for high values of 02 the tt electrons of 
the nitro-group and the benzene ring it is attached to are 
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decoupled, whereas when <j)i is large the tt cloud of the in- 
ner benzene component is increasingly disconnected from 
those of the pyridyl anchors. Both of these disruptions 
of the TT connections between molecular components are 
bound to reduce the effect of NO2 as an electron accep- 
tor, and therefore it is in accordance with expectations 
that a decrease in both, and (/)2, increases EA and EN, 
where the molecular charge AN^I^f is found to follow the 
trend of the electron affinity as the top part of Fig. [3] doc- 
uments at least for all angles above 20°. The downward 
shifts of AN^^f for lower angles (also reflected in the data 
derived for IP) find their origin in strong steric repulsion 
effects for the planar molecule, where we explained the 
reasons in Sec. IIIBI and will illustrate the consequences 
as distortion of the electron density in Sec. IIII El This 
steric repulsion is so destabilizing to nearly planar con- 
formations (see Fig. [TJd) that it puts them out of reach 
for actual experiments and we therefore focus in our main 
analysis and conclusions on angles above 20°-30°. 

We summarize the discussion of Fig. [3] by stating that 
a variation of the angles 0i and 02 can alter the frac- 
tional charge on the molecule by ~ 0.1 electrons. Due to 
the strong Pauli repulsion identified for the pyridyl/ Au 
bond in earlier work^^, AN'^f is negative for a wide 
range of the investigated conformations, which means 
that nitro-bpph becomes positively charged in this junc- 
tion. The conformation dependence of AN^^f is due to 
a decoupling of tt electrons on molecular components for 
increased torsion angles, which is also reflected in values 
for EN and EA calculated for the molecule in isolation. 
We stress that no dependence of AN^^f on (pi exists for 
bpph without a nitro-group (not shown here) which is 
further support for our arguments. 



C. Molecular level positions 

After having addressed the conformation dependence 
of zero bias charge transfer, the next step of our anal- 
ysis is an evaluation of its impact on the alignment of 
molecular levels with respect to Ep. One way of defining 
MOs within the junction in spite of the coupling between 
molecule and electrodes is provided by the NEGF formal- 
isn^, where the Hamiltonian for the electron transport 
is formulated in terms of the localized Siesta basis func- 
tions^ and a subdiagonalization of the molecular AOs 
only can be performed^^. This method gives MO eigenen- 
ergies incorporating all electrostatic interactions between 
the molecule and the leads including charge transfer but 
not direct hybridisation between AOs on separate com- 
ponents of the junction. We show the results of this pro- 
jection scheme for the LUMO and LUMO+1 in Fig. 
in dependence on both and 02- While the level po- 
sitions for the LUMO+1 roughly follow AN^^f and EN 
(as shown in Fig. [3]) in their dependence on molecular 
conformation, the projected values for the LUMO differ 
considerably in their trends. 

In order to investigate this discrepancy further we em- 



ploy another way for predicting the LUMO energy in 
Fig. IDd. Here we make use of the fact that the ener- 
getic position of the lowest lying molecular orbital MOl 
is more than 10 eV below all of the Au valence states 
and is therefore independent of the hybridisation between 
the molecule and the leads^^. This makes it possible to 
read the MOl energy directly from the standard output 
of the DFT calculations for the whole system and add 
the energetic difference between MOl and the LUMO as 
calculated for the molecule in isolation, where both has 
been obtained independently for all values of ^1 and ^2- 
The resulting predictions in Fig. HJ) do not take into ac- 
count any form of direct interaction between the LUMO 
and Au states but they do include all electrostatic effects 
such as charge tranfer and changes in the vacuum level 
alignment due to modified dipole moments (see Sec. IIII El 
for a discussion of the dependence of the molecular dipole 
moment on the torsion angles) as well as variations of the 
energy differences between MOs in the isolated molecule 
due to the distortion. Since these predictions show a 
picture that is in agreement in terms of angle depen- 
dent trends with the LUMO+1 but not the LUMO in 
Fig. H^, we can conclude that the LUMO positions from 
the projection scheme must be governed by the inter- 
action of this particular MO with surface states. This 
finding can be explained by the observation that while 
the molecule can lose electrons due to Pauli repulsion at 
the pyridyl/Au interface, the only way it has of obtain- 
ing fractions of electrons from the Au surfaces is to have 
its LUMO partially occupied (which would automatically 
bring it closer to E^) and that this effect should increase 
with a decrease of both torsion angles (see the discussion 
in Sec. IIII Bp as indeed can be seen in the lower part of 
Fig. Ht. 

Finally, we want to isolate the effect that charge trans- 
fer has on molecular level positioning from other confor- 
mation dependent quantities such as the molecular dipole 
moment and the MO structure of the molecule. For this 
purpose we perform two separate sets of DFT calcula- 
tions on the isolated molecule, where for one of them we 
define a charge neutral setup and for the other we add 
the fractional charges^^ from Fig. [3] to the total num- 
ber of electrons in the unit cell. Assuming that these 
charges are so small that their impact on the position 
of the vacuum potential in the calculation is negligible 
(which appears to be reasonable because the dipole mo- 
ment remains unaltered by the charging), we argue that 
by reading the energy MOl for both sets and forming 
the difference AMOl between them in dependence on 
and (j)2 as plotted in Fig. St, we get a measure of 
the level variation that is caused by charge transfer only. 
This variation is very comparable both in its scope of en- 
ergies (0.2-0.3 eV) and its trends with the torsion angles 
to the conformation dependencies of the LUMO in the 
prediction of Fig.|4]3 and the LUMO+1 in the projection 
scheme of Fig. 2^. From this finding we conclude that 
charge transfer is the dominant source for the conforma- 
tion dependence of the level alignment in the Au/nitro- 
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FIG. 4: (Color online) a) Torsion dependence of the eigenenergies for the LUMO and LUMO+1 as derived from a projection 
scheme for molecular levels within the NEGF formalism'^^ from the electron density of the junction by subdiagonalization of 
the Hamiltonian built from the basis functions localized on the molecule^, b) Variation of LUMO energies with 01 and 02 as 
obtained by adding the energetic difference between the LUMO and the lowest molecular level MOl for the isolated molecule 
to the easily identifiable MOl energy within the composite junction, c) Angle dependent variation of the shift of MOl for the 
isolated molecules due to fractional charging corresponding to the values in Fig. [S] All energy values in this figure are given 
with respect to the Fermi energy of the leads. 



bpph/Au system, because neither the molecular dipole 
moment nor a changed MO structure can have an im- 
pact on the observed variations of AMOl with 0i and 
02 in Fig. [4l3. The charge induced behaviour of AMOl, 
however, mimics the trends in the angle dependent shift 
of the most relevant peaks in the transmission function 
in Fig. [2] which is a central result of our article. 



D. Scissor operator correction for the 
shortcomings of DFT in describing MO energies 

It is well known that the single electron energy gaps 
in DFT are too small when compared to a many body 
description for insulators and semiconductors (including 
finite molecules) due to self-interaction effects and the 
approximate description of electronic correlations^^. In 



addition the gap reduction due to screening effects in 
metallic leads (which are often described by classical im- 
age charge models^^^^^) is also not contained in transmis- 
sion functions calculated from NEGF-DFT at zero bias^. 
Recently, it has been shown that both errors cancel out 
to a high extent at least for the occupied MOs of ben- 
zene molecules adsorbed on metal surfaces^^. Our study 
addresses the relative conformation dependent changes 
of MO energies and their relation to zero bias charge 
transfer^^ within the approximations of a single parti- 
cle NEGF-DFT approach, where so far we have assumed 
that the described shortcomings of our method are not a 
major concern. This assumption was based mostly on the 
argument that the calculated transmission functions can 
be compared with experiments due to the error cancella- 
tion cited above. We stress that the screening effects in 
the leads due to the molecular dipole or interface dipole 
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FIG. 5: (Color online) Shift of the levels clumo and eLUMO-\-i 
in Fig. [2] due to the ad hoc correction with So — ASq intro- 
duced in Ref.^^, where uncorrected values are plotted with 
dashed and corrected values with solid lines. 



moments are properly described within DFT and that it 
is only the polarisation of the metal electrons due to elec- 
tronic excitations of molecular states which are addressed 
in Ref^. 

It is outlined, however, in Ref.^^ that, although there 
will always be some amount of cancellation between self- 
interaction errors and missing polarization effects in a 
semi-local DFT description, the relative size of the two 
contributions might depend significantly on a variety of 
factors such as the structure of the molecule, its orien- 
tation with respect to the surface, the molecule-surface 
distance and the type of substrate. Since some of these 
aspects differ strongly when we compare the systems in 
our work with those of Ref.^^, we want to investigate 
both sources of error explicitly in the following. For this 
purpose and applying the procedure for an ad hoc correc- 
tion in Ref.^^, we introduce a scissor operator — ASq, 
where Hq = —{^lumo-\-EA) emends the energy gap error 
for the isolated molecule and ASq accounts for dynami- 
cal polarisation effects due to the presence of the metal 
surface as calculated from an image charge model (see 
the appendix of Ref.^ for the details of this approach). 
In Fig. Owe plot the uncorrected and corrected clumo 
and €lumo-\-i in their dependence on 0i for 02=0°. The 
following observations can be made: i) The two contri- 
butions Ho and AHq cancel out exactly only for one ge- 
ometry (01=70°). In general Sq varies quite significantly 
with 01 (see also the inset of Fig. [Sb for the 0i depen- 
dence of EA), while ASq = 0.75 ± 0.03 for all struc- 
tures; ii) although the scissor operator correction alters 
the energies of both levels quite strongly in quantitative 
terms, the angle dependent trends discussed in the sec- 
tions above are not only still observable, but their observ- 
ability is even enhanced by the correction, where also the 
corrected levels rise in energy with a decrease of 0i for 
all angles higher than 20-30°. 
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FIG. 6: (Color online) a) Molecular dipole moments fi for 
nitro-bpph in dependence on the torsion angles 01 and 02. b) 
With 02 fixed to 0° , the effect of steric repulsion on the charge 
density is investigated for nearly planar molecules, where the 
density for 0i=9O° is taken as a reference and the sum of the 
charge density is formed over the plane perpendicular to the 
axis defined by the two pyridil-N-atoms. The inset shows the 
structure of the molecule, where the evolution of the total 
electron density along this axis is also indicated in arbitrary 
units. 



E. Molecular dipole moments 

For molecules without a dipole moment in equal dis- 
tance to two identical surfaces, the alignment of molec- 
ular levels is completely defined by a balance of Pauli 
repulsion at the interface, which is the cause for the for- 
mation of interface dipoles in such a setup and the dom- 
inant source for electron subtraction from the molecule, 
and a partial filling of the LUMO, which is the only pos- 
sibility for the molecule to gain electrons. The inter- 
play of these two driving forces for level alignment has 
been investigated for bipyridine and biphenyl-dithiolate 
between Au leads^^^^^. If the junction is asymmetric, the 
picture becomes more complicated. While it has been 
shown that for monolayers of molecules adsorbed on just 
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one surface the level alignment is dominated by inter- 
face dipoles, while molecular dipoles only affect the work 
function^i^, no such rigorous separation of effects can 
be motivated for molecules with an inherent dipole in the 
"sandwhich" setup of the junctions we study in this ar- 
ticle. We face the additional complication that vacuum 
level alignment is also not straightforward in such a junc- 
tion, because two different vacuum levels exist to the left 
and right of the molecule due to the dipole. Throughout 
this article we avoid the latter problem, since we only in- 
vestigate the dependence of level positions on the torsion 
angles 0i and 02, thereby deriving indirect arguments 
for the dominance of charge transfer effects, but do not 
decompose the level alignment explicitly for a given con- 
formation in terms of the relative strength of contributing 
factors. 

In Fig. [6^ we show the molecular dipole moment /i 
as calculated from the electron density and ion distribu- 
tions^^ for all molecular conformations. Strikingly, there 
does not seem to be any dependence of ja on <pi above 30°, 
while only for 02 there are trends which can compare to 
those for charge transfer and level alignment in the pre- 
vious sections. This result seems plausible, because the 
angle 0i determines how much the inner benzene ring 
can attract electrons from the two pyridil groups, where 
there is no reason to believe that this charge transfer phe- 
nomen should exhibit a distinct asymmetry which would 
then induce a change of the dipole moment. The sec- 
ond angle 02, on the other hand, determines how many 
fractions of electrons the nitro-group receives from the 
whole chain of aromatic rings. Since the partial charges 
on NO2 are what causes the molecular dipole moment, 
the relation of /i to 02 becomes evident. This finding 
is additional proof that at least the dependence of level 
alignment on 0i is caused by the change in charge trans- 
fer between the molecule and surfaces, simply because fi 
does not change with 0i . 

Finally, we take a closer look at what causes the dra- 
matic change of /a for 0i below 30° by forming charge den- 
sity differences for these molecular conformations with 
reference to a molecule with perpendicular aromatic rings 
(see Fig. [6)3). It can be seen that the more the conforma- 
tion approaches planarity, the more electrons are pushed 
away from the region where the H and O atoms high- 
lighted in Fig. [TJi meet and partly pushed away from the 
molecule altogether as seen in Fig. [3^1. Interestingly, an- 
other net result of this repulsion effect, is a shift of elec- 
trons towards the right side of the molecule (where there 
is no nitro-group), thereby reducing the molecular dipole 
moment and ultimatively (for 0i=O° and 10°) reversing 



its direction. 



IV. CONCLUSIONS 

In summary, we demonstrated that the shifts of molec- 
ular levels and peaks in the transmission functions in 
the Au/nitro-bpph/Au system induced by a variation of 
molecular conformation can be explained in terms of zero 
bias charge transfer which in turn is governed by the de- 
coupHng of TT electrons between molecular components 
with an increase in torsion angles. The electron accept- 
ing nature of NO2 groups in this junction is countered 
by Pauli repulsion at the interface between the pyridil 
anchors and Au leads, where the latter effect does not 
depend on conformation. 

These findings are relevant for the evaluation of the 
practicability of theoretical device proposals based on ni- 
tro benzene where it is usually taken for granted that 
adding NO 2 to an aromatic ring always induces negative 
partial charges on the molecule and such dependencies on 
molecular conformation and anchor groups are in general 
not discussed or even considered. 

While the torsion angle dependence of the electronic 
coupling in molecules without acceptor groups and rather 
strongly coupled to gold leads has been studied previ- 
ously^^^^^'^^, the effect of molecular conformation on level 
alignment in weaker coupled systems and containing ac- 
ceptor groups has so far not been addressed. The two 
cases differ entirely in the respective mechanisms for the 
conformation dependence. 

Since very recently chemical means have been found to 
fix torsion angles in actual experiments^^^^^, our results 
might open up an avenue for a further exploration of ef- 
fects related to molecular conformation in single molecule 
junctions and even new device schemes for molecular 
switches might result from such experiments. 
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